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Using unitarized Chiral Perturbation Theory methods, we perform a detailed analysis of the 
■KIT scattering poles /o(600) and p(770) behaviour when medium effects such as temperature or 
density drive the system towards Chiral Symmetry Restoration. In the analysis of real poles below 
threshold, we show that it is crucial to extend properly the unitarized amplitudes so that they 
match the perturbative Adler zeros. Our results do not show threshold enhancement effects at finite 
temperature in the /o(600) channel, which remains as a pole of broad nature. We also implement T = 
finite density effects related to chiral symmetry restoration, by varying the pole position with the 
pion decay constant. Although this approach takes into account only a limited class of contributions, 
we reproduce the expected finite density restoration behaviour, which drives the poles towards the 
real axis, producing threshold enhancement and hit bound states. We compare our results with 
several model approaches and discuss the experimental consequences, both in Relativistic Ifeavy 
^ , Ion Collisions and in vr ^ tttt and 7 mv reactions in nuclei. 

^ ' PACS numbers: ll.lO.Wx, 12.39.Fe, 13.75.Lb, 21.65. +f,25.75.-q 
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I. INTRODUCTION 



Oh' 

^ . 

i-S^ . Chiral Symmetry Restoration in a hot and dense environment and its possible experimental signatures are central 
issues in the present program of Relativistic Heavy Ion Collisions and reactions in nuclear matter. In principle, it 
should manifest both in the vanishing of the chiral condensate and in the degeneracy of the possible chiral partners. 
J> ' An important role in this context is played by the cr, the I = J — state corresponding physically to the /o(600) 
I broad resonance listed by the Particle Data Group [l| and observed as a pole in tttt scattering. Since this state shares 
the vacuum quantum numbers, it is a suitable channel to study possible medium effects that may hint towards Chiral 
Symmetry Restoration. An interesting suggestion in this context is the formation of a very narrow state near the 
two-pion threshold as a precursor of the chiral transition. This was first suggested in [2] in a NJL model calculation 
and later on also in the 0(4) model at finite temperature iSj. The presence of such a state was explained in the 
following way. The mass of the a is reduced from its vacuum value towards the mass of the pion by medium effects 
restoring the symmetry. As gets closer to the pion threshold 2m7r (which also changes with medium effects) the 
phase space of the dominant decay cr ^ tttt is squeezed, so that the decay width goes to zero and the physical a pole 
approaches the real axis. As a consequence, there is a threshold enhancement of the a spectral function Such 
an enhancement is produced both by finite temperature and by nuclear density effects [j| and should, in principle, 
^ imply a strong enhancement of the tttt cross section in the I = J = channel 0, Q . The analysis of the tttt scattering 
d [ amplitude in the NJL model shows analogous results at finite density 0] ■ 

More recent analysis have shown some important differences with the works by Hatsuda and collaborators discussed 
above, especially in the finite-T case. A large N analysis in the 0{N) model ^ shows that the physical a pole in 
the second Riemann sheet has still a sizable width at the temperature for which the real part of the pole has already 
reached 2mTr. Further increasing of medium effects does reduce the imaginary part to zero. At that moment, two 
poles coexist on the real axis: one of them is the original second-sheet a pole and the other one comes from its 
unphysical partner in the upper half plane. From there on, one of the poles moves onto the first Riemann sheet, 
becoming a stable bound state. In fact, there had been earlier proposals about the formation of a tttt bound state 
in nuclear matter as a consequence of the increasing in-medium strengt h [p . This behaviour is argued in to be 
universal, provided there is a stable state above the critical point. In pTlT^ a self-energy study of the a pole at 
finite T in the 0(4) model is also presented, but the results are slightly different from Q. These authors find a real 
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pole in the second Riemann sheet for all temperatures, coexisting with the physical <j pole. The real pole moves with 
T towards threshold, causing threshold enhancement when it arrives 2m^. However, the enhancement is smeared out 
if the pion thermal width is taken into account [lH. From that temperature onwards, that pole crosses to the first 
sheet, becoming a finite temperature tttt bound state, while the physical a pole remains on the second sheet, reducing 
its real part towards 2mj but keeping a nonzero imaginary part and its broad nature. Finally, it is worth mentioning 
also the analysis in [isl . [3| of nuclear effects in tttt scattering using effective lagrangians in a chiral unitary approach, 
where a weaker threshold effect than in [F] is obtained, and the importance of p-wave pion renormalization in tttt 
scattering is stressed. 

Experimentally, there is evidence for a / = J = tttt threshold effect induced by the nuclear medium in ttA tttt A' 
[isl [l6j and tttt A' [TtI reactions in nuclei. The CHAOS collaboration [l^ has reported a much stronger effect 

than the Crystal-Ball one There is then some controversy both in the results and in their interpretation. Thus, 
while in [l5| a peak in the two-pion invariant mass distribution is seen at threshold for increasing nuclear density, 
the results being compatible both with chiral-symmetry restoration and ordinary nuclear effects in the a channel, the 
effect in (l6| a nd [l3| is weaker and seems to be better described by models in which the a appears as a broad tttt 
resonancefisj . No such threshold nuclear effect is seen in the isospin channels 1 = 2 and 1=1 . 

In this work we will use unitarized Chiral Perturbation Theory (ChPT) to investigate the behavior of the pion 
scattering poles as the system approaches chiral symmetry restoration, with the sole physical input of chiral symmetry 
and unitarity. Although ChPT provides a perturbative expansion in meson energies and temperatures and therefore it 
is not capable to predict a true chiral phase transition, its predictions are based only on the chiral symmetry breaking 
pattern and are therefore model independent. In fact, the quark condensate calculated within ChPT shows a restoring 
behaviour, predicting a critical temperature around ~ 200-250 MeV when it is extrapolated from the low-T region 
On the other hand, the requirement of exact unitarity within the ChPT framework has led to a fruitful 
unitarization programme of the meson- meson scattering amplitudes [2ll . [2^ . [23l . [24| , one of its main results being the 
generation of all the low-lying meson resonances with their masses and widths in agreement with the values quoted 
by the PDG [if. Recently, the pion scattering amplitudes have been calculated and unitarized at finite temperature 
[25I [26|. which has allowed to determine the low-T evolution of the /o(600) and p{770) poles. Here, we will study 
those poles, attending to their relation with Chiral Symmetry Restoration. The specific new aspects that will be 
studied here and were not considered before are: i) Precursors of Chiral Symmetry Restoration such as threshold 
enhancement, emphasizing the modifications of the original arguments when a not-narrow state is present (section 
IIII|) , ii) A detailed analysis of poles in the real axis below the tttt threshold, both in the second (virtual states) and in 
the first (bound states) Riemann sheets (section lIVp . This will require a generalization of the unitarized amplitudes 
in order to account properly for the presence of Adler zeros and remove unphysical poles, which is also a new result in 
vacuum, iii) The nature of the /o (600) and p states as revealed by their thermal and in- medium behaviour (sections 
IVl and IVip and iv) The analysis of nuclear density chiral restoring effects (section IVip . 



II. FORMALISM 



Chiral Perturbation Theory is the most general low-energy framework compatible with the SUL{Nf) x SUu{Nf) 
SUv{Nf) Spontaneous Chiral Symmetry breaking pattern of QCD with Nf = 2, 3 light quark flavors. We will restrict 
to the two flavor case here. ChPT is built as an expansion in p/A^, p denoting generically any pion momenta or mass 
and A^ ~ 4:71/.^ ~ 1 GeV is the typical chiral scale. Considering the heat bath temperature T as a formally 0{p) 
quantity, ChPT can also be used to study a meson gas at temperatures below Tc [isl . fioj . 

The original calculation of pion-pion elastic scattering amplitudes in ChPT to one loop was given in [27| and 
its finite-T extension is discussed in [2^ . Consider a partial wave t^'^{s;T) with definite isospin / and angular 
momentum J, where s is the center of mass energy squared of the pion pair. The generic ChPT structure is then 
t^'^{s; T) = t2'^ {s)+tl"' {s; T) + . . . , where tk denotes the 0(p/ K^)^ contribution. The lowest order ^2(3) is temperature- 
independent and gives Weinberg's low-energy theorem [28] from the lowest order chiral lagrangian £2 (the non-linear 
sigma model). It has the following form: 



ti\s)=A"{s-sl,') (1) 

Here, sq is the so called Adler zero (the point where the amplitude vanishes) to second order. For the relevant 
partial waves in pion scattering, the values of the A and sq constants are given in Table [J Note that to this order, 
the amplitudes depend only on the pion mass m^r — 140 MeV and the pion decay constant ~ 93 MeV. The term 
includes contributions both from tree level £4 and one-loop £2 graphs. The tree level contributions are T-independent, 
renormalize the amplitudes to this order and depend on certain low-energy constants. We follow the convention in 
[2?! and write the amplitude in terms of the li constants, which are independent of the renormalization scale. The 
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TABLE L Values of the constants for the partial waves at lowest order in ChPT. 

one-loop pion scattering amplitude depends only on li and I2, but if it is written in terms of the physical and 
as we do here, then it also depends on ^3 and I4. The numerical values of the low-energy constants are fixed 
phenomenologically. We will use the same values as in [2^, namely li = —0.3, I2 = 5.6, I3 = 3.4 and I4 — 4.3, which 
fit the mass and width of the p(770) (see below) when the partial waves are unitarized. 

The one-loop part in contains in particular the imaginary part required by unitarity. At T = 0, exact unitarity 
for partial waves reads lmt^'^{s) = cro(s)|i^'^(s)p or, for the inverse amplitude: 



Im 



= -ao{s) ^ t^'^is) = 



Re [l/t^^is)] - iao 



(2) 



above the two-pion threshold, i.e, for s > 4m^ where cro(s) = ^1 — 4m2 / s is the two-pion phase space. The ChPT 
expansion cannot reproduce unitarity exactly, but only pcrturbatively, i.e, Imt4 — (To(s)|t2(s)P and so on. Typically, 
the unitarity violations are more severe as energy increases, especially if there are physical resonances with the same 
quantum numbers as the channel under consideration. For these reasons, several unitarization approaches have been 
proposed, which amount to approximate Re[l/t] in ([2|) in different ways. A simple and successful one within the 
ChPT context is the so called Inverse Amplitude Method (lAM) [2l|, (2^ which is built by demanding exact unitarity 
and matching with the ChPT expansion up to fourth order when expanded at low energies, i.e, t^^'^'^ = t2 + + . . . , 
where the different chiral orders can be traced by counting powers of 1//^. We recall that the perturbative expansion 
is meant to converge at T = near threshold so that it is in that region where the formal expansion in inverse 
powers of has a physical meaning. With these two requirements, the lAM amplitude reads for a given partial wave 



t'^^' = [t2?/ih-U). 

At finite T, the partial waves calculated in p| 



satisfy the following thermal perturbative unitarity relation 



Imi4(s;T) = aT(s)|i2(s)r 



(3) 



where: 



(4) 



and nsix) — {exp{x/T) — 1)^^ is the Bose-Einstein distribution function. 

The function axis) is the thermal phase space [295, which is increased with respect to the T = one by the 
difference [1 -|- nsiEi)] [1 -I- n_B(£^2)] — nB{Ei)nB[E2) — 1 + nB{Ei) +nB{E2), where Ei^2 are the energies of the two 
colliding pions. The first term in the difference corresponds to the enhancement of the available phase space due to 
the increase of two-pion states by the scattering of two pions in the thermal bath (these are the stimulated emission 
processes discussed in [29]). The second term reduces the phase space by absorption due to collisions of the incoming 
pions with the thermal bath ones. In the center of mass frame, Ei = E2 = \fsl1 and the thermal phase space reduces 
to (j4|). Therefore, one can extend the usual T = exact unitarity arguments, replacing (Tq — > ctt and the partial 
waves by the finite-T ones, provided that only intermediate two-pion states are relevant in the thermal bath. This is 
reasonable if we remain in a dilute gas regime, which is consistent at low and moderate temperatures. Thus, from the 
perturbative amplitude in the center of mass frame, following the same steps as for T = case, the unitarized lAM 
thermal partial waves for every J J channel are given by [2^ : 



(5) 
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Therefore, the above thermal unitarized amphtude satisfies the exact unitarity relation ([2]) with (Tq replaced by ctt 
and the partial waves replaced by their finite-T counterparts obtained in (25j . As a consequence, it has the same 
analytic structure as the T — one, i.e., it is analytical on the complex plane off the real axis, it has the right unitarity 
cut starting at the two-pion threshold s — 4m^ and a left cut for s < coming from diagrams in the t — u channels 
[2^ . From these analytical properties, one can immediately define the amplitude on the second Riemann sheet which 
continuously connects with the imaginary part of the first Riemann sheet amplitude t^-^*^ in ([5|) across the right cut, 
i.e, Imi^^(s — ie) — Imt^(s + ie) = aT{s)\t^\'^ ~ — Imi^(s — ie) for s > 4rnJ. One has: 



' l-2iaT{s)ti{s-T) 

The poles of t^^ in the lower complex half plane correspond to physical resonances. We follow the standard 
convention and denote the pole position by Spoie = {Mp — iTp/2Y [i| chosen so that Mp and Tp would correspond 
to the mass and width of a narrow Breit-Wigner resonance. Since we are working by definition in the center of 
mass frame of the incoming pions, the masses and widths correspond to resonances at rest. At T 0, and with the 
low-energy constants chosen above, t^^'^^ has a pole at Mp ~ 441 MeV and Tp ~ 464 MeV, corresponding to the 
/o(600). Note that this state is clearly not a simple narrow Breit-Wigner resonance, since its width is of the same size 
as its mass. In the I = J ~ 1 channel, the pole is at Mp ~ 756 MeV, Tp ~ 151 MeV, corresponding to the p(770) [l|. 



III. THERMAL EVOLUTION OF THE POLES 



At r 7^ 0, the lAM poles are shown in Figure [TJ The p pole follows a finite-T trajectory compatible with dilepton 
data in Relativistic Heavy Ion Collisions, namely, a broadening behaviour with only a small mass decreasing. This 
point will be commented in more detail in section [V] The thermal widening in the p channel can be explained in 
part due to the increase of thermal phase space given by cft- In fact, it can be shown (26j that if the change in 
the mass and in the effective pirn coupling constant are neglected (see below) the width would increase simply as 
Tt/To = CTT (M^) /(To {Mj) = [1 + 2nB{Mp/2)], which is rather accurate up to T ~ 100 MeV. For higher T, there is 
a further increase of the width due to the increase of the effective vertex [2611. 
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FIG. 1: Evolution with temperature of the unitarized /o(600) and p complex poles in the second Riemann sheet. The points 
are obtained by varying the temperature in 20 MeV intervals. 



The phase space argument cannot be directly applied to the /o(600) pole, since this is a broad resonance. Still, 
there is a broadening effect for moderate temperatures, as it can be seen in Figure [1] However, this behaviour changes 
qualitatively when the temperature is further increased, due to the effect of Chiral Symmetry Restoration, which is 
notorious in this channel. The system tends to restore chiral symmetry by reducing drastically the mass and in this 
sense this state tends to become degenerate with the chiral pion triplet (see section |V|. In principle, the argument 
proposed in 0, Q would suggest that as Mp 2m^, the width should go to zero by phase space reduction. However, 
it is very important to note that this is a narrow resonance argument and it does not necessarily hold for a wide 
state, as it is the case here. In fact, observe that our pole remains as a wide state even at temperatures close to the 
phase transition, where its mass has already reached the two-pion threshold. In fact, due to the initial phase space 
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widening, the width close to the transition is of the same order as for T = 0. This is in agreement with the behaviour 
found in Q befow the transition and it is one of our relevant results. Thus, we do not expect this pole to cause any 
sizable effect in the real axis near threshold, like an enhancement of the cross section, which could be regarded as 
a precursor of chiral symmetry restoration. This is indeed the case, as observed in Figure [21 Note that almost no 
variation with temperature is seen at threshold, i.e, the scattering length remains almost T-indepcndcnt. This result 
had also been obtained in the finitc-T scattering length calculation in [30] and a qualitative explanation is provided 
below. To fully confirm the absence of chiral restoration precursors at finite T, we need to perform an analysis of the 
real poles, which deserve a separate treatment for the reasons explained in section HVl 
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FIG. 2: Squared modulus of the t unitarized partial wave for different temperatures. 



Let us now provide a qualitative description of the above results. Consider the propagator and spectral function of 
the resonance under consideration at rest: 



s~AP{s) + tM{s)T{s) ^'^^ 
Pis) = -2lmD ^ '-^^i^ , (8) 

Recall that ([8]) reduces to p = 2TrS{s — M^) when F ^ 0^, as it corresponds to a free particle of mass M and 
positive energy. Here we are denoting by M^(s) and [— M(s)F(s)] the real and imaginary parts of the resonance 
self-energy respectively, with both M and F real functions of s in the general case. In the complex s plane, the pole 
of the propagator is the solution of Spoie = M'^(spoie) — 'iM{spoie)^{spoie)- Thus, with our convention for the poles, 
we have: 



p2 ^2p2 
M^{.Spole)^M^^-^ ; T'ispole)^J^^^^ (9) 

Consider now the contribution to the tttt scattering amplitude of the exchange of a resonance R with propagator 
(O and the same quantum numbers as the corresponding IJ channel: 



s — M''(s) + iM (sjl (s) 

where Ma is the effective ttttR vertex squared. The parametrization PH)) of the partial waves can be made unitary 
by demanding ([2]) in (jlOp . which gives F(s) = cro(s)a(s) for s € M and s > Arn^. This parametrization is therefore 
compatible with chiral symmetry and unitarity, although a correct description of the 00 channel in the real axis would 
also need the addition of background contributions [Slj, which does not change our qualitative arguments here. As far 
as we do not make any claim about the relative size of F and M nor about their behaviour with s, the parametrization 
pO)l can be viewed as a generalized Breit-Wigner form on the real axis above threshold. 

It is instructive at this point to remind the narrow resonance limit behaviour. If Fp ^ Mp, the pole is close to 
the real axis, so that the amplitude on the real axis is peaked around s = AP ~ Mp and the width Fp ~ F(Af^) = 
a{IvP)aQ{AP)6{AP — Anri^). If M ^ 2m7r, driven for instance by chiral restoration in the 00 channel, the spectral 
function IS]) and Imi are then enhanced for s near threshold as: 
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a^j s — 4to^ 

, , 2m7r6' is - Ami) , ^ 

Imi(s) (12) 

This is the typical threshold enhancement produced by phase space squeezing [s^, '3] which consequently implies an 
enhancement in the tttt cross section, as explained before. On the other hand, if medium effects do not change M and 
a. much, as it happens for the p, the only temperature modification is the replacement of the phase space ctq ^ fT, 
which leads to r^/ro = (tt/'^Oi announced above. 

Let us now show how the previous arguments change when the resonant state is not narrow, i.e, when the pole 
position values Tp and Mp are of the same order. For that purpose, we will consider, as a working example, a BW-like 
form (fTOj) for the amplitude. 

First, in the broad resonance case, the dependence of the functions M, a, F with s may be important. For instance, 
if M is roughly constant and large compared to F, a formal s-expansion of (fTO]) gives t ^ a/M . Comparing this to the 
lowest ChPT order t2 in ([T]) (required by chiral symmetry) would give a ^ AAI{s — sq) near threshold. If sq 4m^ 
as in the 00 channel, that behaviour for a is compatible with the enhancement as M — > 2TO7r, since near threshold 
the denominator of (jlOp is dominated by its imaginary part F ~ ctq a-s compared with the real part (s — M^) ^ CTq. 
However, this need not be the case if M as given by (O is not large compared to 2m7r (for the /o(600), M ~ 373.3 
MeV already at T = and it gets further reduced at finite T). In that case, t(s) in (fTO| near threshold can be made 
compatible with t2 by taking a ^ — A{s — sq){s — M'^) / M , so that t{s) ~ {— M a) / {s — A'P) ~ A{s — sq) near threshold, 
even if AI —> 2m7r at a certain temperature, since now F ~ ctq so that the ctq in the real part of the denominator 
dominates. The result is now that t{s) near threshold is roughly independent of M and therefore of T, with no 
threshold enhancement. This is consistent with what we see in Figure [21 which means that the leading order ChPT 
expression (T-independent) basically dominates at all temperatures near threshold. 

Second, the phase space vanishing of the width near threshold for narrow resonances relies heavily on the fact that 
the spectral function is peaked around s ^ Af^ or, in other words, that the pole is reasonably close to the real axis. 
However, for the /o(600), the spectral function is broadly distributed so that the phase space axis) is not directly 
evaluated at s = M^. An explicit way to see this is to consider the decay rate of the R ^ tttt process. In the narrow 
resonance limit, R can be considered a particle with four-momentum = so that the differential decay rate in 
the rest frame is given by: 



where 



^ n (^^^(2^)''^^'^ [p^pi- (14) 



is the two-particle differential phase space, P = {Mii,0), Ef = + and therefore s ~ [pi +^2)^ — P^ = M^. 
Replacing the i? — > tttt T-matrix element in (fT3)) by Ma, which is assumed to depend only on the energy of the R state, 
and integrating the phase space J d$i2 — (To{s)9{s — 4m^)/(87r) gives for the decay rate F^ — cro(M^)a(M^)6'(s — 
4m^)/(327r) — F(A/^)/(327r) ~ Fp/(327r), yielding the standard relation between the width and the decay rate, which 
is also suppressed by phase space when the mass of the state approaches threshold. 

If the spectral function of the resonance is not well approximated by a 5 function, the differential decay width (jl3[) 
is generahzed by replacing — > J{dui/2TT)p{ui)9{ui) so that instead of Tjy one has the number of decays per unit 
volume and per unit of four- momenta [s^l , integrated over the energy lu distributed according to the spectral function 
of the resonant state. Accordingly, with the spectral function ([5]) and replacing M,r,ao,a by their T-dependent 
counterparts, the quantity 

f°° ds MT{s)aT{s)aT{s) 2Mt{s)Tt{s) 

i'T — I 7: Tr—F= 9 ^ (15) 

74^2 27r 2^5 [s - M^{s)f + [MTis)rTis)f 
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should behave quahtatively as a generahzed decay rate and therefore we expect it to follow a similar pattern as the 
imaginary part of the pole position Tp{T). As a consistency check, in the narrow resonance limit Fx — aT<JT0{M!^ — 
4m^). In the general case though, when Mr — > 2rn^, there will be a non-vanishing contribution to the integral due 
to the broad spectral function. 

Obviously, we cannot draw any conclusion about the evolution of Ft without specifying the T and s dependence 
of ax, Mt,Tt- In order to check that this is a consistent description, let us choose for simplicity s-independent 
AIt and Ft, related to the pole position values by with the values of Fp and Mp obtained in our full lAM 
calculation. We also take ax = F^, which is consistent with the unitarity requirement away from threshold (the 
spectral function in the integrand is suppressed for s near threshold in the broad case, according to our previous 
arguments) . We emphasize that this is a rough approximation, since the actual profile of the partial wave does 
not resemble a BW-like form, even with large F [3ll]. Nevertheless, it provides us with a simple model to explain a 
thermal resonance behaviour similar to the one we obtain. With that simplification. Ft rises with T until T ~ 100 
MeV {Fioo/Fo — 1-1) due to phase space increasing in the numerator of p3|) and then starts decreasing as the real 
part of the pole approaches threshold, reaching a value -F20o/^o — 0.7. Comparing with the evolution of Fp in Figure 
III where Fp(100 MeV)/Fp(0) ~ 1.2 and Fp(200 MeV)/Fp(0) ~ 0.8, we see that a broad resonant thermal state with a 
finite generalized decay rate near the two-pion threshold behaves in a qualitatively similar way as our thermal fo (600) 
state. 

From the above discussion we conclude that for a broad state like the /o(600) at finite temperature there is always 
a balance between thermal phase space increasing, which has nothing to do with chiral symmetry restoration, and 
chiral restoring mass and width decreasing. Depending on how broad the resonance is when chiral restoring effects 
start being important, the pole may move towards the real axis more or less rapidly. There is no way to know a priori 
whether it can be close enough to be considered a narrow resonance before the chiral transition, so that observable 
effects such as threshold enhancement may take place. In our approach, this does not happen at finite T (the situation 
may change at finite density as we will see in section rvT| i.e., the /n (600) is still broad at the transition, in agreement 
with Q but not with previous approaches [1, 0] ■ We also agree with [O, as far as the complex broad / = J = 
pole is concerned, and in the next section we will check for the presence of real poles. 



IV. REAL AXIS POLES AND ABLER ZEROS: 



The analysis of poles near the real axis below threshold is relevant for the chiral restoring behaviour, as discussed 
above. As we will show here, the unitarized lAM amplitude below threshold has to be extended in order to account 
properly for the Adler zeros and remove spurious poles. For clarity, we will discuss the T = case first. 



A. T = extended amplitude 

Consider the lAM unitarized amplitude (O. Near threshold, t^^^^ ^ ^2+^4 and ^4 is perturbatively small compared 
to ^2- On the other hand, the perturbative amplitude must have an Adler zero, i.e, a point where it vanishes, as a 
consequence of Weinberg's low-energy theorem |28j , which predicts that the pion scattering amplitude vanishes in the 
chiral limit at s = 0. Both t2 and ^2 + ^4 have this property and is natural to expect that t^^'^^ has a zero not far 
from the perturbative one. The zero of ^2 is given by sq in ^ and we will call si the zero of ^2 + ^4- Note that the 
lAM expression (O vanishes exactly at sq, although the order of the zero is not necessarily the same: if ^4(50) 7^ 0, 
^lAM ^ gj^")2^ while t2 has a zero of order one. In addition, note that if t2 and ^4 have the same sign and ^4(50) 7^ 0, 
there is a point s > sq where the denominator of ([5]) vanishes. If that point exists for a given channel, it produces a 
spurious pole in the first Riemann sheet, since it would predict a not-observed T — tttt bound state. That pole does 
exist in the / = J = channel at s ~ (89 MeV)^ and at s ~ (197.4 MeV)^ for / = 2, J 0. In the / = J = 1 channel 
there is no such problem, since the scattering length has to vanish by parity so that sq = si = 4rn^ and at that point 
both t2 and ^4 vanish, so that ([5]) is not divergent. The position of sq, si and s (all of them below threshold) for the 
00 channel are shown in Figure [31 The contribution of the possible amplitude zeros was already mentioned in early 



lAM papers [24], although it was not taken into account since its effect is negligible in the physical region, as we 
will show below. The existence of the non-physical pole has also been noticed in [s^, where it is also linked to the 
presence of the Adler zero. 

Therefore, the lAM definition ([5]) is affected with two related problems below threshold: the wrong order of the 
Adler zero and the presence of spurious poles. We will see now a simple way to cure both without spoiling the 
essential properties of chiral symmetry and unitarity of the lAM. For that purpose, we note that if we extend the 
inverse amplitude as: 
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FIG. 3: Unitarized amplitudes at T = in the 7 = J = channel, where so and s\ are, respectively, the Adler zeros of second 
and fourth order and s is the point where the lAM diverges. 



^71^-^7^+.9(-) (16) 

with g an analytic function off the real axis, real for s £ M, the unitarity condition remains unaltered and so does 
the analytic structure of the amplitude. The choice of g is in principle arbitrary. Therefore, since having a double 
zero in the amplitude means having a double pole in the inverse amplitude, if we choose 5fc(s) — ~Rl{s — sq)^ with 
R the residue of 1/t at the double pole so, we get an amplitude with an Adler zero of order one at s = sq. In fact, 
note that, taking into account ([1]), the expansion of the inverse amplitude near sq reads: 



t^A^is) A2(s-so)2 A2(.s-So) ^ ' 

where the dots stand for terms which do not diverge as s — > sq. 

Therefore, R = —t4{so)/A'^ and the extended amplitude with this choice of g becomes simply: 



fiAMbf^ = [^2(^)]^ ns) 

t2{s)-Uis)+Uiso) 

The expression ([T5| has a simple Adler zero at s = sq as long as ^4(50) 7^ 0. But in addition, if ^2(5) — ^4(s) is an 
increasing function of s from s = to threshold, the JAM pole disappears, since in that case, sq is the only point 
where the denominator of ()18p vanishes. This is what happens in the 00 channel, as shown in Figured 

We can improve further the unitarized amplitude, by demanding also that it matches the ChPT series near threshold. 
Note that this condition does not hold for (US)), which can easily be seen by reexpanding it in powers of f^'^. Since 
^4('So) = CI/tT**), we get t{s) ~ ^2(3) + ^4(5) — ^4(50) instead of the ^2 + ^4 of ChPT. That is, it only matches correctly 
the first 0{f^^) order ^2- If we want the amplitude to match ^2 + ^4) then, as we will show below, it is sufficient to 
demand that both have the Adler zero at the same point, i.e, at s = si. Therefore, in order that 1/t has a simple 
pole at that point an according to (|17p we must subtract now the sq pole and add the si one, i.e, we choose: 



5c(s) 



1 



s - So 



s - Sl 



(19) 



where c is an undetermined constant that we will fix by demanding the perturbative matching with the ChPT series 
to fourth order. Adding HI]) to we get: 



^lAMc^. ^ A^{s-So)\s-Si) 

C(s - So)2 - (S - Sl) [t4{s) - t4{so) - {S - So)t4{so)] 



(20) 



9 



Let us expand the previous expression in powers of /~^. Recall that A = 0{f~'^), = 0[f~^) and Si = so+C(/~^), 
so that the leading order is t^^^'^ — {A^ /c){s — Sq) and therefore c = A + 0{f^^). Expanding now the expression 
^2(51) + ti{si) = around sq we find: 



Si = .so-i4Gso)M + 0(/r) (21) 

Using this in the expansion of to fourth order, we have t^^^^''{s) = t2{s)+t4(s) - {s- so){c- A+t'4^(so))+0{f^^). 
Therefore, by taking: 



c^A-t'^iso) (22) 

in (j20p , the unitarized amplitude matches the chiral expansion up to fourth order and has the Adler zero at the same 
position and with the same order as the perturbative amplitude. Note also that with c in ((22|) . (|20|) reduces to ([18]) 
when si — + so, as it should. In fact, since the difference between si and so is perturbatively small, we expect this 
extended amplitude to behave very similarly as the simplest version psp . and in particular to remove the spurious 
pole while being numerically closer to the original lAM amplitude. That is indeed the case, as it can clearly be seen 
in Figure [3l The same situation takes place in the / = 2, J = channel. The extended amplitudes are expected to 
differ little from i^^^^ away from sq or si, since the g^.c functions vanish for s ^ Sq, Si. For that reason, in Figure[3] 
we have only shown the region around the spurious pole and the Adler zero. Away from that region, either in the real 
axis or in the complex plane, there is practically no difference between the lAM or its extended versions discussed 
here. 

It is interesting to note that if we had not redefined the amplitude, there would be also a non-physical pole in the 
second Riemann sheet just below the real axis and below threshold. The lAM in the second Riemann sheet for T ~ 
reads, from ^ and (O: 



^lAMJI ^ (23) 

t2 — ^4 — LCrto 



where we have defined d[s) — icro(s — iO"'") = Am^ /s — 1 for < s < 4m^. Thus, the denominator of (|23p is positive 
near threshold for an attractive channel (dominated by t2 > 0) like the 00 one and diverges to minus infinity as 
5^0+, so that it must have an odd number of zeros. Since the denominator is not zero at s = sqi that means that 
there would be at least one such real pole for the lAM in the second sheet. Consider however the extended amplitude, 
for instance ^TE\\ . Now, 



^IAMb,II _ [2 /2^^ 

which is not necessarily singular below threshold since both the denominator and the numerator vanish at s = sq. A 
similar argument holds for ^^^^.a/cJ/^ whose denominator is now positive both at threshold (provided again that the 
0{f^) terms dominate there) and at s = 0. The conclusion is then that for the extended amplitudes, there are no 
real poles in the second Riemann sheet below threshold, or an even number of them. We have checked that there are 
no such real poles at T = 0. It is unclear whether the real second sheet pole found in [ll|, [l^] (for the a self-energy) 
is also spurious, in the sense discussed here. 

Finally, we have checked that the p(770) and /o(600) poles remain at the same place at T = with the second 
Riemann sheet extensions of the extended amplitudes, as expected from our previous arguments. We also note that 
the above extended amplitudes can also be obtained by using dispersion relations [s^ . as it was also done in early 
derivations of the lAM f22'|, which provides a formal justification of the results shown here. 



B. T 7^ corrections 



Once we have extended the T — amplitude to ensure the correct behaviour below threshold, we obtain the 
extended T amplitudes by including the thermal corrections in i4(s; T). Therefore, the values of si and s become 
temperature dependent. With this extended amplitude, we do not find any spurious pole in the first Riemann sheet, 
as it happened for T = 0, up to temperatures close to the phase transition, nor any additional pole in the second sheet 
apart from the standard /o(600) pole. The results are shown in Figure 2] with the extended amplitude i^'^*^^(s; T). 
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We remark that this is consistent with the fact that the /o(600) pole shown in Figure [T] remains far from the real axis 
near the phase transition. Our result in this respect is different from that in [Til. Il2j|. where an additional pole in the 
second sheet is found for all temperatures, although we agree with [ll[ in that finite-T threshold effect are not seen. 
We also agree with where the additional pole does not appear for low and moderate T. 
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FIG. 4; Extended unitarized amplitude at finite temperature, in the first and second Riemann sheets. No real pole is observed 
for temperatures up to T ~ 200 MeV. 



An important remark is that temperature effects make i4(s;T) grow with respect to t2^s) so that some of our 
previous T = conclusions could change. For instance, regarding the number of poles, the denominator of (jlSp or 
(PU)) could become now negative at threshold, so that we would have an odd number of real poles in the second 
Riemann sheet and also in the first one, which would correspond to tttt bound states. Although this is not the case 
at finite temperature, we will see in section IVTl that it is possible that other chirally restoring effects, such as nuclear 
density, may cause the appearance of new poles in both sheets, as those obtained for instance in 8] and consistently 
with threshold enhancement nuclear effects. 



V. THE NATURE OF THE THERMAL RESONANCES 



The analysis of tttt scattering poles at finite temperature or density may also be useful to understand the nature 
of the resonant states in terms of quark and gluon degrees of freedom. This is particularly relevant for the /o(600), 
for which there is a long-standing controversy about its nature. In 0(4) or linear-a models, the ct is introduced as 
an explicit degree of freedom in the lagrangian, and by definition it transforms as a pure gq, isospin singlet state 
[sH - On the other hand, the /o(600) listed by the PDG [l| is described by several approaches, including unitarized 
ChPT [13, H^, US] > where resonant states are dynamically generated, i.e, not included as explicit degrees of freedom 
from the start, nor making any requirement about their nature. When strangeness is included, unitarized ChPT 
generates precisely the scalar resonances corresponding to a scalar nonet pattern (/o(600), /o(980), ao(980) and the 
controversial k) as well as the standard vector meson octet [H, [l^, HI] . Thus, these results support the existence 
of such a scalar nonet, in agreement with models where those states are included explicitly. However, there is an 
important discrepancy about their qq nature. In fact, a recent analysis of the large-A^c behaviour of the low-lying 
resonances generated by unitarized ChPT [s^ shows that poles of the members of the scalar nonet, such as the 
/o(600), do not behave as a pure qq state unlike for instance the p or the other members of the vector meson octet. 
The observed /o(600) state main content could then be of mm-qq states with the same quantum numbers, such as 
tetraquark (or two- meson) or glueball states [s^, [s^l ■ 

Therefore, when medium effects are included in models where the a is explicitly introduced [1, 0, S HHi i 
expected that m„ ^ (a), the order parameter, as temperature or density is increased [12]. In that framework, it is not 
surprising that near the phase transition, the a becomes degenerate with the pion, both in mass and width, since in 
the chiral limit the only source for the a mass is its vacuum expectation value so that the complex a pole approaches 
the real axis. However, as we have commented above, the observed /o(600) might not have such a simple structure, 
and this should show up in the pole behaviour as the system approaches the chiral phase transition, although bearing 
in mind the limitations of our approach as T approaches Tc- 

Our unitarized ChPT results show that the /o(600) exhibits a clear tendency to reduce the mass of the pole as 
the temperature increases. For T greater than about 100 MeV, this effect makes also the width decrease, overcoming 
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the width increase due to thermal phase space. This could be interpreted as a tendency of this state to become the 
singlet state partner of the pion in the limit of chiral symmetry restoration. However, the size of the thermal width 
when the chiral transition has already been reached indicates that this state is not degenerated with the pion. Thus, 
although a state with the numbers of the qq isospin singlet survives near the critical point, it does not seem to behave 
like a qq state at finite temperature. 




T (MeV) T (MeV) 

FIG. 5: Comparison between the /o(600) pole relative mass and width with the quark condensate, for the physical pion mass 
(left) and near the chiral limit, for m,r = 10 MeV (right). 



The situation is summarized in FigureEJ where we compare the pole mass Mp{T) with (g(7)i/2. The reason for doing 
so is that {a) in the 0(4) model corresponds to /,r to leading order. On the other hand, from the Gell-Mann-Oakes 
Renner (GOR) relation = —mq{qq) /m'^ [111 where niq = = md, we expect /j- to scale with T like {qq)^^^ if 
the T-dependence of m-^ is ignored. This is a good approximation in one-loop ChPT 18] whereas for two-loops there 
are pion mass corrections to the GOR relation when it is expressed in terms of T-dependent quantities 39] . We have 
taken the quark condensate from the virial expansion [20| , where the pressure and the quark condensate are expressed 
in terms of T = scattering amplitudes. The results shown in Figure [5] are obtained with the 0{p'^) pion scattering 
amplitudes, giving a critical temperature Tc ~ 230 MeV, similarly to the ChPT perturbative three loop calculation in 
19] for the pion gas. The differences introduced by considering unitarized amplitudes are small in the virial approach 

m]- 

In Figure [5] we also show the results for = 10 MeV, i.e, close to the chiral limit. Here, a technical remark 
is in order. As commented above, we are expressing our amplitudes in terms of the scale-independent low-energy 
constants li, with the same convention as in [27|. Namely, they are related to the one-loop renormalized constants as 
= l^ [h + log /(327r2) where 7i are numerical factors and ^ is the renormalization scale. Since the /f (^) 

are mass- independent, when we change the pion mass from mi to TO2, we should also change h h -I- log {ni\/m'^ in 
the amplitudes when calculating the poles. The same observation applies to the numerical value of /tt we are using, 
which is also mass-dependent since to one loop, = / [l -|- /{IGn^ p)^ [23], where / ~ 87.6 MeV is the pion 
decay constant in the chiral limit. 

For the physical pion mass, the mass curve departs from the condensate one as it approaches 2m7r. This is consistent 
with the idea that when the qq mean value is negligible, the non-qq component of the /o(600) state still survives in 
the thermal bath as a short-lived state. The departure from the cr-like evolution is even greater if effects beyond 
one-loop are considered. For instance, m^(T) calculated to first order in the pion density [40] with the lAM T = 
partial waves, decreases considerably near the chiral transition so that it gets even further from Mp{T), while f-n{T) 
at two-loops starts increasing from T > 150 MeV [39|. 

The above conclusions are confirmed when reducing the pion mass, so that explicit chiral symmetry breaking terms 
are negligible, chiral symmetry restoration takes place at a lower temperature (for ttItt = 10 MeV we get Tc ~ 182 
MeV, also close to the chiral limit three- loop result in [l^) and the cr-like behaviour for the pole should be more 
visible. The results shown in Figure [5] for to^t = 10 MeV are qualitatively not very different from the physical mass 
case. The values of the mass and width of the pole at T = are now Mp ~ 406.2 MeV and Fp ~ 522.7 MeV. 
Therefore, the pion mass reduction produces a larger effect in the width (due to the increase of phase space) than 
in the mass, consistently with the idea that the main contribution to the /o(600) mass comes from the spontaneous 
chiral symmetry breaking scale and not from the explicit part. The larger width makes the interpretation of the 
/o(600) in terms of a real particle still more doubtful. At finite T, the width increases again first by thermal phase 
space and decreases near the transition along with the pole mass, which decreases monotonically. We see that the 
mass curve is now closer to the condensate than for the physical pion mass, which confirms our previous comment 
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about restoration in the chiral limit. However, the width remains again very large at the critical point and in fact the 
mass is even far from 2to^, revealing once more the possible presence of non-gq components. The GOR deviations 
and pion mass T-dependence are negligible near the chiral limit [s^, 113] . 

These results, combined with the absence of real poles for the extended amplitudes discussed in section ITVl which 
we have also checked for TOtt — 10 MeV, lead us to conclude that in our finite T approach there is no /o(600) — tt 
degeneracy at the critical point. 

Finally, we comment on the implications of our results for the finite temperature nature of the p resonance. As 
discussed in section Hill our main result is that near the critical point the p exhibits a substantial thermal broadening, 
which increases with temperature. The mass reduction is small, of only about 50 MeV at T = 200 MeV. This 
behaviour is compatible with other model calculations and its implications for the pion electromagnetic form 
factor [4l|, which is the quantity entering directly into the dilepton rate, are consistent with recent experimental data 
on the CERN-SPS dilepton spectrum [il, Q and on the p shape from the tt+tt" spectrum in RHIC-STAR 45]. 
On the other hand, the original Brown-Rho scaling hypothesis [4a] predicts that the p mass should decrease with 
medium effects roughly as the condensate, i.e, Mt/Mq ~ {Q'i)t / {qq)o, at least close to the transition point. This 
hypothesis has been supported by approaches based upon the hidden local symmetry framework 47], where the p is 
introduced as a gauge field of the local chiral symmetry, and implies that at the critical temperature the p should 
become degenerate with the pion in the chiral limit. The BR dropping-mass scenario predictions [i^ are compatible 
with the CERES dilepton data [43] but they seem to be in conflict with the recent NA60 dimuon data [3, E^] ■ The 
latter has been nevertheless questioned in [47, 50], where, among other arguments like the violation of Vector Meson 
Dominance, it is claimed that the BR Mp dropping takes really place only near Tc, whereas the mass changes very 
little until about T ~ 125 MeV. 
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FIG. 6: Comparison between the p pole and the chiral condensate. Here, we include a simple fit of the form suggested by 
dropping-mass scenarios, where To ~ 310 MeV is the temperature at which the extrapolated p mass vanishes. 



We present in Figure [6] our results for the p mass compared to the quark condensate, where we have extrapolated 
our JAM poles shown in Figure [T] to the temperature To ~ 310 MeV where the pole mass would vanish. Note that 
Tq is notably far from the chiral transition point estimated by the vanishing of the chiral condensate at Tc ~ 230 
MeV, calculated in the virial approach as commented previously. We have also represented a fit of our results with 
a curve of the type Mp{T)/Mp{Q) = [1 — (T/Tq)"]" with n integer, as suggested in [49j. We get a good agreement 
for n = 6 and a ~ 0.8. Therefore, although our results are compatible with a dropping mass behaviour with only 
sizable mass decrease very near To [13, [13] , the value of such "critical" temperature is rather high and the mass curve 
does not follow the quark condensate. On the other hand, as commented above, our main effect is the p broadening, 
whereas in [i^] the width does not change much except near the critical point where it also drops. In addition to 
the thermal phase space increasing, this is related to the behaviour of the puir effective vertex which also drops in 
(47| . whereas in our case grows with T (26| . making the width grow even more. In conclusion, unitarized ChPT is 
not fully incompatible with the BR scenario in what concerns the qualitative behaviour of the p mass, but the nature 
of the p state when the system approaches the chiral transition is very different in both approaches: we find a wide 
state of relatively little change in mass, while in BR-like approaches the p tends to a massless narrow state. In other 
words, the thermal behaviour of the p in our approach is dominated by non-chiral restoring effects, unlike in the BR 
scenario. We also point out that, due to the vector nature of the p resonance, it is not obvious that one can relate 
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its chiral restoring behaviour with its qq content, as we did in the scalar /o(600) case which has the vacuum quantum 
numbers. 



In the previous sections, we have analyzed the temperature behaviour of the tttt scattering poles in the 00 and 
11 channels. One important conclusion of that analysis has been the notorious influence of thermal phase medium 
effects, which in the /o(600) case compete with chiral restoration with the result that no observable effect is produced 
when the real part of the pole approaches the two-pion threshold. This thermal phase space effect is the result of the 
scattering with pions in the thermal bath, as explained in section ITll 

In this section we will consider the influence of other chiral restoring effects such as nuclear density. For that 
purpose, we will follow a simplified treatment, taking T = and encoding those effects in an effective decreasing 
The GOR relation is valid to linear order in nuclear density at T = [5l| as long as the /jr associated to the time 
component of the axial current is used (it is generally different from the space component in the medium "^^l) and 
therefore, we expect to scale like the quark condensate in the nuclear medium. In addition, the density dependence 
of the pion mass is rather weak and can be neglected in a first approximation so that to leading order [51[ : 



where p is the nuclear density, cj-^m — 45 MeV |53| is the pion-nucleon sigma term and po — 0.17 fm^"^ is the normal 
or saturation nuclear matter density. 

In this approach the interaction between pions and nuclear matter is considered in an effective mean-field way. This 
is the same philosophy followed in earlier linear-cr model analysis 0, [H] which introduce the linear density dependence 
only through the function — (a) (p) / (a) {0) ~ /^(p)//^(0) . The same procedure of simulating medium effects 
for pion scattering only in the reduction of /^r is followed in [1, Q . More elaborated finite-density analysis of the 
pion scattering amplitude (3] conclude that to linear order in density the amplitude corrections roughly amount to 
consider the in-medium pion-propagator and pion decay constant. We do not consider corrections in the pion self- 
energy, but even when those are taken into account |3] the result is still a significant threshold enhancement with 
no in- medium broadening of the /o(600). The approach followed in [141] includes additional in- medium corrections 
and the dropping behaviour of is automatically incorporated [11]. In that approach the pion propagator is also 
renormalized in the medium, other effects being subdominant, and the results obtained for threshold enhancement are 
similar to other approaches. All these works support the conclusion that for low nuclear density and T — there are 
no intermediate states producing a significant broadening as with the thermal space broadening discussed in section 
im which is proportional to Bose-Einstein pion distributions and is therefore absent at T = 0. For instance, in [3], 
the renormalized in-medium pion propagators induce effectively new intermediate states which can be interpreted 
as coming from collisions with matter particles, but these new channels do not make the /o(600) broader [ssj . 
Nevertheless, we stress that we do not mean to describe the full tttt scattering at finite density only with this approach. 
Our main purpose is to see whether pure chiral symmetry restoration terms as those given by (|25p can reproduce 
typical threshold enhancement effects in our unitarized ChPT scheme. 

Although we expect the reduction of to be the dominant effect for precursors of chiral symmetry restoration 
in the / = J = channel, this is not so clear for the I = J = 1 channel, for which no observable threshold effects 
are seen in tttt scatterin g in nuclei [13]. The BR scaling suggests that the mass of the p resonance should also drop 
with increasing density |46l . [isl . [s^ and near the critical density it should follow the condensate curve (pS)) . Since 
the p resonance is a narrow state, this implies that the width also vanishes by phase space reduction, so that the 
pole is driven by chiral restoration towards threshold, as the /o(600) pole. The crucial point, as far as finite density 
effects are concerned, is whether the threshold enhancement is visible for densities p<po, which are those accessible 
experimentally [Tsl. [la. [l7| . 

The evolution of the poles with decreasing is shown in Figure [T] In both channels, we have used the extended 
amplitude (PO)) . since the distance of the pole to the real axis is reduced for increasing density. A salient feature is 
that in both channels the poles move towards the origin. This can be qualitatively understood by noting that the 
lAM (|23)) scales with as t = i2/{f^i2 — i4 — ai"^), where the functions with tilde are independent of and similarly 
for the extended amplitudes. The domain of convergence of the chiral expansion (/|i2 3> ?4,?| around threshold is 
therefore reduced when decreasing and thus the value at which the denominator vanishes may be closer to the real 
axis. 



VI. CHIRAL RESTORING MEDIUM EFFECTS 
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FIG. 7: Evolution of the complex poles with decreasing /,r at T = 0, corresponding to chiral-restoring finite density effects 
through (|25p . Attached to every point is the value of /,r in MeV units. All the points correspond to poles in the second Riemann 
sheet, except the black triangles which are in the first sheet (bound states). 



The /o(600) pole mass reaches threshold at about ~ 55 MeV, or density p ~ 1.86po from (pS)) . Since the phase 
space effect is absent, the resonance becomes narrow so that the mass reduction impUes also a reduction in the width, 
which rapidly drops to zero when the real part reaches threshold. The effects on threshold enhancement are notorious 
as we show in Figure[S]for = 70, 60, 55 MeV which according to (P5|) correspond to densities p/po = 1.2, 1.7 and 1.9 
respectively. Therefore, we reproduce qualitatively the threshold enhancement at finite density obtained in previous 
works 0, H, [E 0) 9 ■ The quantitative values are not far either: in the enhancement of jt^^p at threshold for 
$ = 0.7 (/tt ~ 78 MeV) is about 9 times its vacuum value, while we get a factor of 6.5 for /^r — 70 MeV. In 0], 
Mo- ~ 2m-^ at a density p ^ 2po, in good agreement with our valueSj^hile in Q the complex pole reaches threshold 
at about f-„ ^ 65 MeV. The results for the pole mass and width in [55| are also similar to ours. 
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FIG. 8: Squared modulus of the t'^^ and t^^ unitarized partial waves for different values of 



It is even more interesting what happens when we let the density increase further. At a given point (/tt — 52 
MeV) the pole width vanishes and there are two real poles in the second sheet, which correspond to virtual states. 
Remember that, according to our discussion in section HVl the number of real poles in the second Riemann sheet has 
to be zero or an even number once we have extended properly the amplitude to remove spurious poles and provided 
the amplitude remains positive at threshold. When the density is increased a bit further, one of those poles jumps to 
the first sheet, becoming a tttt bound sate. This happens in the 00 channel at f-^ — 45 MeV, as indicated in Figure [T] 
For that value, the extended amplitude at threshold becomes negative because the 0{f^) t2 term no longer dominates 
and therefore the same argument of section IIVI leads to the conclusion that there is an odd number of poles in the 
first and in the second sheets, one pole in this case. This pole behaviour is in agreement with that found in [8|, where 
the "pole doubling" in the second sheet is found at about /tt ~ 60 MeV and the movement of one of the poles to the 
first sheet is at /^r ~ 50 MeV. 

The evolution of the /o(600) pole with density in this approach is now closer to a tr-like qq term, as shown in Figure 
ini The mass follows now quite closely the curve of /^r for increasing density. We are plotting the second sheet pole. 
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which we observe that approaches the pion mass and in fact it does it more rapidly than the /,r curve, approaching the 
/J one. Our conclusion is then that at sufficiently high densities a virtual qq-\ike state, or any other state behaving 
in the same way near chiral restoration, degenerates with the pion and coexists with a tttt bound state. Nevertheless, 
we should bear in mind that our approach relies on the ChPT power counting and therefore it losses validity as is 
reduced. 
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FIG. 9: Scaling with of the pole masses in the / = J = (left) and I — J — 1 (right) channels. The points are the mass 
values of the second-sheet poles in Figure [T] the dashed line is the linear density dependence f-rrip)/ f-wiO) and the full line is 
the quadratic /^(p)//^(0) which scales approximately like the quark condensate, c.f., ea. l|25p . 



As for the I = J = 1 channel, first we see that the pole is still rather away from the real axis when the /o(600) 
pole arrives threshold, which would produce a less important threshold enhancement for this channel for the relevant 
density range, as it can be observed in Figure[8l consistently with the previously mentioned experimental information. 
When the density is further increased, there also appears a bound state joint by a virtual state in this channel, 
although the arguments used in section HVl for the /o(600) channel do not apply here since the perturbative amplitude 
vanishes at threshold. This simultaneous "softening" of the /o(600) and p poles for decreasing f^^ has also been noted 
in @. 

We also see that if only these finite density effects are considered, the p mass follows quite closely a BR scaling 
pattern, as represented in Figure [5] The mass scales with /^r, i.e, with {qq)^/"^ and gradually moves to the {qq) curve. 
This is in agreement with recent analysis in this context [5y| which show that the vector meson masses scale like 
{qqV^'^ for low densities and like the condensate near the transition point. In any case, we insist that this would 
reflect only the finite density effects associated with chiral restoration, which may not be the only relevant medium 
effects. 



VII. CONCLUSIONS 

Using methods based on the unitarization of the elastic tttt scattering in one-loop Chiral Perturbation Theory, we 
have studied the behaviour of pion scattering poles with the presence of medium effects such as temperature and 
density, as the system approaches the chiral phase transition. The unitarization framework we have employed is the 
Inverse Amplitude Method. For the temperature evolution we have used thermal amplitudes, whereas we mimic chiral 
restoration effects associated to nuclear density by studying the variation with the pion decay constant in the linear 
density GOR approximation at T = 0. 

The I — J = pole in tttt scattering near the chiral phase transition remains a broad state with a pole position 
Mp ~ 2to^ but nonvanishing Tp. This implies that only with finite temperature effects is very unlikely that there 
could be observable precursors of chiral symmetry restoration like an enhancement of the tttt scattering cross section, 
as claimed by various authors. As we have shown in detail, the reason for that is that thermal phase space increases 
the width considerably before chiral restoring is effective, so that near the critical point the resonance spectral function 
is widely spread. In fact, we have shown that our results can be qualitatively understood in terms of the generalized 
decay rate of a wide resonance into two pions. 

At finite temperature we do not find additional poles in the second Riemann sheet near the real axis (virtual states), 
which is consistent with the /o(600) remaining a broad resonance. For the analysis of real poles, it is crucial to extend 
the unitarized amplitude in order to account properly for the Adler zeros. In this way, we have shown that we can 
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get rid of non-physical poles appearing below the tttt threshold in the 00 and 20 channels. The extended amplitudes 
arc almost indistinguishable from the lAM ones in the physical region, i.e, above threshold. This part of our work 
is interesting also for formal reasons, since one docs not lose the powerful unitarity and chiral symmetry constraints 
when demanding also a behaviour below threshold consistent with the perturbative Adler zeros. 

We have provided an interpretation of our results at finite T in terms of the non-degeneracy of "chiral partners" 
near the critical temperature. The fact that the /o(600) width remains large at the critical point, even near the chiral 
limit, supports this fact and is also consistent with a sizable presence of non qq components in the thermal /o(600) 
state. This is confirmed also by the differences between the mass and the quark condensate curves as T approaches 
the critical point. Our results for the p state in the thermal case predict a broadening behaviour with slight mass 
reduction and are compatible with several broadening-p scenarios. They are in conflict however with dropping mass 
descriptions, with which we only agree in the dropping behaviour of the mass, but with a vanishing temperature larger 
than the critical temperature predicted by the quark condensate. 

When the pion decay constant is reduced, to reproduce linear density chiral restoring effects at T = 0, we obtain a 
sharp chiral restoration pattern in the I = J = channel, showing threshold enhancement at densities similar to those 
obtained in different models and qualitatively consistent with experiments, in spite of the fact that our finite density 
approach includes only a limited type of contributions. When the /o(600) pole reaches threshold, two poles appear 
on the second Riemann sheet, one of them remains a virtual state that tends to degenerate with the pion, while the 
other becomes a two-pion bound state. The nature of the /o(600) within this approach resembles much more the 
typical qq one of linear a models. As for the p, the finite density effects associated to chiral symmetry restoration 
produce a Brown-Rho like scaling of the pole mass, which moves towards threshold first as the square root of the 
condensate for moderate densities and as the condensate for higher ones. The p gets narrower with only these effects 
and a bound- virtual state pair also appears at high densities. 

In summary, the presence of chiral symmetry precursors related to real poles in the J = J = channel is disfavored 
in our approach when finite temperature effects arc included, while pure density effects favor these precursors and 
Brown-Rho scaling. Therefore, in a Relativistic Heavy Ion Collision, where both temperature and density effects are 
relevant, our approach does not predict any measurable effect related to the pole evolution other than the p widening, 
observable in dilepton data. On the other hand, our results for density-like effects are compatible with in-medium tttt 
production experiments, where threshold enhancement has been observed, despite the limitations of our approach. 

As a future extension of this work, a finite temperature and density analysis of the meson-meson scattering poles 
in the SU (3) case may also reveal important information about Chiral Symmetry Restoration, where the presence of 
heavier states becomes important. 
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